Soliton motion in a parametrically ac-driven damped Toda lattice 
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We demonstrate that a staggered parametric ac driving term can support stable progressive 
motion of a soliton in a Toda lattice with friction, while an unstaggered drivng force cannot. A 
physical context of the model is that of a chain of anharmonically coupled particles adsorbed on a 
solid surface of a finite size. The ac driving force models a standing acoustic wave excited on the 
surface. Simulations demonstrate that the state left behind the moving soliton, with the particles 
shifted from their equilibrium positions, gradually relaxes back to the equilibrium state that existed 
before the passage of the soliton. Perturbation theory predicts that the ac-driven soliton exists if 
the amplitude of the drive exceeds a certain threshold. The analytical prediction for the threshold is 
in reasonable agreement with that found numerically. Collisions between two counter propagating 
solitons were also simulated, demonstrating that the collisions are, essentially fully elastic. 



The important role of collective nonlinear excitations 
in the form of solitons in several condensed-matter 
physics contexts is widely recognized [jlj. Particularly, 
discrete nonlinear models have received an increased 
amount of attention lately . Often actual observation 
of nontrivial dynamical effects in realistic systems are im- 
peded by friction. Compensation of friction and support 
a soliton therefore usually require an externally applied 
driving force ||. In the simplest case, progressive motion 
of a soliton in the presence of dissipation is supported by 
a dc driving force Q , which in fact is the only possibility 
in a homogeneous continuum medium. A less obvious 
but physically important option is to drive a soliton by 
an ac force (with no dc component) either in a continuum 
system subjected to a periodic spatial modulation ||, or 
in a discrete system in the form of a chain of nonlinearly 
interacting particles HQ. 

A simple example of the latter system is the ac-driven 
damped Toda lattice (TL) 



_i) - exp {x n+ i - x n ) = 
-ax n + (— l) n e sin(wt) , 



(1) 



where x n is the displacement of the n-th particle, e and 
uj are the amplitude and frequency of the driving force 
and a is the friction coefficient. This model describes 
a chain of nonlinearly coupled particles in a dissipative 
environment driven by an externally applied electric ac 
field. The alternating sign in front of the driving term, 
(—1)" (a discrete pattern of this type is frequently called 
staggered j^]), implies alternating electrical charges of the 
particles in the chain; without the multiplier (—1)™, the 
driving term can be trivially eliminated from Eq. (1). 
In Refs. [||, it was demonstrated analytically and nu- 
merically that the model (Q) indeed supports progressive 
motion of a soliton in the presence of finite friction (this 
was accomplished for this model in the so-called dual 
form). Ac-driven motion was predicted and found at the 
resonant soliton velocities 



V = oj/tt {2v + 1) , v = 0, ±1, ±2, . 



(2) 



The objective of the present work is to study, moving 
solitons in a similar but essentially different model, viz., 
a parametrically ac-driven damped TL, with the driving 
taken both in the staggered form, 

x n + exp (x„ - x n -i) - exp (x n+ i - x n ) = 

— ax n + (—l) n e sin (ut) ■ x n , (3) 

and in the unstaggered form, differing from Eq. (0) by 
the absence of the multiplier (—1)™. The model (|J) de- 
scribes a physical system in the form, for example, of a 
chain of particles adsorbed on a solid surface and inter- 
acting through an anharmonic repulsive potential. Ad- 
ditionally, the particles interact with a uniform exter- 
nal field aligned perpendicular to the adsorbing surface. 
(The flield may simply be gravitation, or a dc electric 
field). In the latter case, all the particles are assumed to 
have the same charge (and the interaction between them 
within the chain is generated by the charges through 
Coulomb repulsion). The effective drive is induced by 
a standing acoustic (elastic) wave excited on the surface, 
the end particles in the chain being fixed at the surface 
edges. If the length of the surface (i.e. the length of the 
chain) is L, and the number of the particles is N + 1 , 
in the equilibrium state the repelling particles are sep- 
arated by the distance L/N. The finite size of the sur- 
face selects the wavelengths of the standing elastic waves 
A = 2L/M, where M is an arbitrary integer (the number 
of the standing half- waves). Commensurability between 
the adsorbed chain and the standing wave takes place 
provided M — pN, where p — 1,2,3, ... . If the com- 
mensurability index p is odd, starting from p = 1, the 
distance between neighboring particles in the adsorbed 
chain is equal to an odd number of the standing half- 
waves, hence the ac driving forces acting upon the neigh- 
boring particles are it out-of-phase, yielding the staggered 
from of the driving term in (|^). Oppositely, when p is 
even, the phase shift between the neighboring particles is 
a multiple of 2ir corresponding to the unstaggered driv- 
ing force. Finally, it is necessary to specify the boundary 
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conditions for the standing elastic wave. If the edges of 
the surface are fixed, the particles in the adsorbed chain 
are located at the nodes (zero-amplitude points) of the 
standing wave. It is easy to see that in this case the 
external field oriented perpendicular to the surface gives 
rise to the direct driving- force term in Eq. (1). On the 
other hand, if the edges are free, the particles are lo- 
cated at the standing-wave maximum-amplitude points, 
which gives rise to the parametric drive in Eq. (for 
sufficiently small x n ). 

In the following we shall first numerically demostrated 
that the model (0) is able to sustain propagating soliton- 
like solutions which velocity locked to the frequency of 
the drive. Further, we demonstrate that only in the pres- 
ence of the staggered driving force is the soliton like so- 
lutions supported. The locking of the velocity to the 
driving frequency and the fact that the soliton like solu- 
tions needs a staggered drive to exist is then explained 
analytically and the threshold of the driving amplitude 
is found. Finally the relaxation dynamics of the particles 
in the wake of the soliton is studied analytically. 

The system of equations (1) was solved numerically 
for a chain consisting of 500 particles. The natural ini- 
tial condition to apply is the exact soliton solution of the 
unperturbed TL 



x„ (t) = - In 
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(l + £-2{n-n -Vt)) 



(4) 



where the real parameter £ (— 1 < £ < 1) determines the 
properties of the soliton (width and velocity), and n is 
an arbitrary phase constant. The velocity V is given as 



V=(£-l/0/hi(r 2 ) 



(5) 



Launching the soliton (Q) into the system (@) gives rise to 
scenarios such as those depicted in Fig. 1(a) and 1(b). In 
the case of Fig. 1(a) the friction is very small (a = 0.01) 
and the initial kink transforms into a propagating front 
and in the wake of the front the particles suffer relatively 
strong oscillations and only slowly relax to their equilib- 
rium position. In the case of stronger friction (a = 0.18) 
Fig. 1(b) shows that the front is still created but the os- 
cillations of the wake particles are strongly overdamped 
and the particles therefore relax rapidly to their initial 
positions. In the latter case the net shift generated by 
the passing front is zero, which in some sense makes the 
driven front "more solitary" than the soliton of the un- 
perturbed TL. 

The situation is quite different when an unstaggered 
driving force is applied. Figure 2 shows a representative 
example of the motion of a driven front in this case. The 
front is seen to form in much the same fashion as in the 
model with the staggered drive. The front is however 
decreasing in amplitude and eventually ceases to exist. 
However, it can exist for rather long times, up to 100 
periods (2-k/cS) of the driving force. 

Since the front in the staggered model behaves simi- 
larly to a soliton, an interesting experiment is to launch 



counter propagating fronts into the same system and ob- 
serve their behavior. An example of such an experiment 
is shown in Fig. 3 and it is clearly seen that the fronts 
survive the collision. A closer study of trajectories of the 
two fronts shows that the velocity is unchanged by the 
collision and that no phase shift occurs. 

We now investigate the observed phenomenon analyt- 
ically and specify the region of existence in parameter 
space. We will show that the velocity of the propagating 
front is controlled by the frequency of the driving force. 
As in ||, the analytical consideration will be based on 
perturbation theory, assuming a and e to be sufficiently 
small. In the zeroth-order approximation, the exact soli- 
ton solution to the unperturbed TL equation is given in 
Eq. (|J). Passage of the soliton through a given particle 
gives rise to a displacement 



Ax n = 



c n (t)dt = sgn£ • In (<f 



(6) 



If the drive is taken in the staggered form, as in Eq. (||) , 
the analysis similar to that of Ref. || shows that the 
ac-driven motion of the soliton is expected at the same 
resonant velocities (^). The unstaggered drive, however, 
gives rise to the resonant velocities 



V = cj/2-kv, v = ±1,±2, 



(7) 



The goal of the analytical approach is to predict the 
minimum (threshold) value e^ nr of the ac-drive ampli- 
tude e that is sufficient to compensate the friction and 
allow progressive motion of the soliton. Each value of 
the integer v in (g) and (0) should thus give rise to a 
corresponding function e+^i- (a, co). Steady motion of the 
soliton through the damped driven lattice is possible if 
the momentum lost by each particle under the action of 
the friction force, OLXji • IS balanced by the momentum 
input from the driving force, (— l) n esin (uit) ■ x n . The 
momentum loss of the particle can be calculated as 



(A J P) k 



+ 00 



a / x n (t)dt = — sgn£ • aln (£ 



(8) 



where Eq. (H) was used. The momentum input is 



(AP) 



input 



(-l) n e / sin. (cot) ■ X n (t)dt 



= (-iyw 



+ OC 



cos (cot) ■ X n (t)dt, (9) 



where integration by parts was performed and the con- 
tribution from t — > ±oo was neglected (since it is zero on 
avarage ). The integral in (g) can be calculated explic- 
itly, inserting the expression for X n (t) following from the 
unperturbed soliton solution (Q). The result depends on 
the arbitrary phase constant no; the threshold is deter- 
mined by equating the maximum possible absolute value 
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of the momentum input (g) to the absolute value of the 
momentum loss (|^). After a simple calculation, this leads 
to the equality 



2tt 



aln (£ 2 ) = — e thr • 



sinh 



(10) 



Replacing the combination (£ — 1/£) _1 In (£~ 2 ) by what 
follows from Eq. (^|), and using the resonant relations 
(2) or (Q), one immediately concludes that in the case 
(7), corresponding to the unstaggered drive, the expres- 
sion multiplying e* n j. in Eq. (HG) vanishes. This implies 
that the threshold (^hr) becomes infinite, simply indi- 
cating that the progressive motion of the soliton cannot 
be supported by the wnstaggered ac drive in the presence 
of friction, which agrees with our numerical observation. 
However, substitution of Eq. (||) for the staggered driv- 
ing leads to the meaningful result 



e thr = 2^ ln (6> : 



sinh 



Co - 1/& 



«/(«), (ii) 



where the parameter £o is determined, as a function of v 
and u>, by a transcendental equation following from Eqs. 
@ and |) 



7T (21/ + !)(&-!/&)= win (£o" 2 ) 



(12) 



Note that the proportionality of to a in Eq. jTl| ) 
is a trivial feature of the first-order approximation of the 
perturbation theory while the dependence e+hr^) 
is a nontrivial issue. Equations ( pi] ) and ([l2]) give this 
dependence in an implicit form. To obtain an explicit 
dependence, one needs to solve Eq. ( |l2] ) for £o, which 
cannot be done analytically in an exact form, but an ap- 
proximate solution is available, 



Co" 1 



2lu 



tt(2is + 1) 



ln 



2lu 



tt(2is + 1) 



(13) 



provided that the driving frequency is sufficiently large, 
so that 
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tt(2z/+ 1) 



> In 
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7t(2z/+1) 



(14) 



In this case, the expression ( |TT[ ) takes an approximate 
explicit form, 



e thr 



auj 
= — ln 

7T 



2lo 



7T(2z/ + r 



x sinh 



^tt 2 (2i/ + l)/ln. 
2 V 77 \.7r(2i/ + l) 



2w 



(15) 



Numerical simulations show that the propagating soli- 
ton can be created and sustained at the resonant veloc- 
ities given by Eq.(||) when the the soliton parameter £ 
is chosen in accordance with Eq.(n2). Further we find 



that the soliton velocity locks to the driving frequency 
even when the initial velocity deviates slightly from the 
resonate velocity. That is if the initial velocity fulfills 
(lo — 5u>) /n < V < (lo + Slo) / ir, where Slo ~ 3, the soliton 
will be created and sustained at the resonant velocity. 

The threshold characteristic of the driven motion was 
measured numerically in the following way. At a fixed 
value of the driving amplitude e, the simulations were 
run at different values of lo, gradually increasing the fric- 
tion coefficient until reaching a maximum value a max , 
beyond which the drive could no longer support the sta- 
ble moving soliton. In terms of Eq. ([Ll]) , 



"max = e/f (to) 



(16) 



i.e., this procedure implies a way to numerically measure 
the function /(lo). 

In Fig. 4, the numerically obtained dependence 
"max (w) is plotted, again with the analytically predicted 
dependence (0), at the fixed value of the amplitude 
e = 0.2. The analytical dependence is shown in two 
forms: the simplified Eq.(15) (dashed line), and the more 
accurate form (dashed-dotted line) obtained by numeri- 
cally solving the transcendental equation ( |l2|) and sub- 
stituting the result into Eq. ( |TT| ) . The parts of both an- 
alytical curves beneath the turning points are irrelevant: 
a direct inspection shows that they correspond either to 
unphysical solutions with £ 2 > 1, or to very broad soli- 
tons, for which, in fact, the model becomes cwerdamped, 
and the perturbation theory does not apply. A conclu- 
sion suggested by Fig. 4 is that, although the accuracy 
provided by the perturbation theory is not very high, 
it captures systematics reasonably. It is also noteworthy 
that the simplified approximation ( |l5| ) works much worse 
than the full perturbative approximation based on Eq. 
( |l2| ) this is not surprising: even for the largest frequency, 
to = 22, for which the result of the direct simulations is 
given in Fig. 1, the condition (|l4|), necessary for appli- 
cability of the simplified approximation, takes the form 
2.63 > 0.97. 

Finally, the slow relaxation of the oscillations behind 
the passing soliton can be easily analyzed. Indeed, as- 
suming the oscillation amplitude to be small enough, one 
can linearize Eq. (ph, yielding 



\%7l— 1 2x n -(- Xn+l) 

( — l) n e sin (Lot) x n — ax n 



(17) 



A solution to Eq. ([17]) can be sought as a sum of a 
rapidly oscillating staggered component and a slowly re- 
laxing unstaggered one 



x n (t) = (-l) n x(t) + X(t). 



(18) 



Substituting this into ( |l7j ) and collecting the staggered 
and unstaggered terms, we obtain two simple equations: 



x + ax + Ax — eX sin (Lot) , 
X + aX = e sin (cot) ■ x. 



(19) 
(20) 
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Using the fact that, in the present approximation, the 
friction is much weaker than the inertia, and following 
the assumption according to which the unstaggered part 
of the solution is slowly relaxing (in comparison with the 
rapid oscillations of sin(wt)), one can readily obtain the 
following approximate solution to Eq. (^o|): 



it) 



-4 



sin (ut) ■ X(t). 



(21) 



Next, substituting ( pi| ) into (g(j) and replacing the 
rapidly oscillating term with its avaraged value, we ob- 
tain an effective equation governing the slow relaxation 
of the state left behind the passing soliton: 



X + aX + -<? - 4) 



X = 0. 



(22) 



Depending on a relation between the small parameters a 

and e 2 (u 2 — 4) , this equation may describe both oscil- 
latory and aperiodic relaxation to the final state X = 0, 
which is the equilibrium state that existed before the 
passage of the soliton. Furthermore, a numerical study 
of the relaxations shows that the frequency specified by 
Eq.(p2[) is the relaxation frequency observed in the full 
simulation. 

In summary we have demonstrated, analytically and 
numerically, that a staggered parametric ac driving term 
can support stable progressive motion of solitons in a 
damped Toda lattice, while an unstaggered drive can- 
not. Also, the threshold condition for the existence of the 
ac-driven soliton predicted by our perturbation theory 
is in reasonable agreement with that found numerically, 
though there remains a clear difference. The simulations 
demonstrate that the state left behind the moving soli- 
ton gradually relaxes back to the equilibrium state that 
existed before the passage of the soliton, which we eas- 
ily explained analytically. Finally, we demonstrated that 
collisions between two solitons moving with opposite ve- 
locities are nearly elastic. 
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FIGURE CAPTIONS 



Fig. 1: Typical examples of a moving front with a stag- 
gered driving force, (a) is with small friction (a = 0.01) 
and (b) is with relatively strong friction (a — 0.18). In 
both cases e = 0.2, and ui = 12.4). 

Fig. 2: Representative example of a moving front 
with an unstaggered driving force. The parameters are 
a = 0.02, e = 0.2, and u = 12.4). 

Fig. 3. Typical example of the quasielastic collision 
between two ac-driven solitons. The values of the pa- 
rameters are e = 0.2, a — 0.1, and ui — 12.4. 

Fig. 4. The dependence between the driving frequency 
lo and the maximum friction coefficient a at which a sta- 
ble moving soliton can be supported by an ac drive with 
fixed amplitude e = 0.2. The dots are results of direct 
numerical simulations of Eq. (|^). The dashed-dotted 
and dashed curves represent, respectively, the analytical 
results provided by the full perturbative approximation 
based on Eqs.fll]) and (|i~2|), and by the simplified ap- 
proximation ([15). The parts of the curves beneath the 
turning points are irrelevant (see text). 
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